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The complete set of analytic solutions of the geodesic equation in a Schwarzschild-(anti-)de Sitter 
space-time is presented. The solutions are derived from the Jacobi inversion problem restricted 
to the set of zeros of the theta function, called the theta divisor. In its final form the solutions 
can be expressed in terms of derivatives of Kleinian sigma functions. The different types of the 
resulting orbits are characterized in terms of the conserved energy and angular momentum as well 
as the cosmological constant. Using the analytical solution, the question whether the cosmological 
constant could be a cause of the Pioneer Anomaly is addressed. The periastron shift and its post- 
Schwarzschild limit is derived. The developed method can also be applied to the geodesic equation 
in higher dimensional Schwarzschild space-times. 


PACS numbers: 02.30.Hq, 04.20.-q 

I. INTRODUCTION AND MOTIVATION 

All solar system observations and almost all other 
observations related to gravity are perfectly described 
within Einstein’s General Relativity. This includes light 
deflection, the perihelion shift of planets, the gravita¬ 
tional time-delay (Shapiro effect) the Lense-Thirring 
and the Schiff effect related to the gravitomagnetic field, 
as well as strong field effects governing the dynamics of 
binary systems and, in particular, binary pulsars 
However, there are two phenomena which do not fit into 
this scheme and still represent a mystery; that is dark 
matter and dark energy. Dark matter has been intro¬ 
duced to explain the galactic rotation curves, gravita¬ 
tional leasing, or particular structures in the cosmic mi¬ 
crowave background. Dark energy is needed to describe 
the accelerated expansion of the universe. All related ob¬ 
servations like the fluctuations in the cosmic microwave 
background, structure formation, and SN la are con¬ 
sistently described by an additional energy-momentum 
component which appears in the Einstein field equation 
as a an additional cosmological term 

Rfiu ~ (f) 

where A is the cosmological constant which, using the 
independent observations mentioned above, has a value 
of |A| < 10“^^ m“^. 

As a consequence, it is necessary in principle to de¬ 
scribe all observations related to gravity within a frame¬ 
work including the cosmological constant. However, due 
to the smallness of the cosmological constant it seems 
unlikely that this quantity will have a large effect on 
smaller, that is, on solar system scales. In fact, it has 
been shown within an approximation scheme based on 
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the frame given by the Schwarzschild-de Sitter space- 
time that the cosmological constant plays no role in all 
the solar system observations and also not in strong field 
effects [3,13 . Also within a rotating version of this so¬ 
lution, the Kerr-de Sitter solution, no observable effects 
arise Q- Nevertheless, there has been some discussion on 
whether the Pioneer anomaly, the unexplained accelera¬ 
tion of the Pioneer 10 and 11 spacecraft toward the inner 
solar system of apioneer = (8.47 ± 1.33) x 10“^^ m/s^ 
0 which is of the order of cH where H is the Hubble 
constant, may be related to the cosmological expansion 
and, thus, to the cosmological constant. The same or¬ 
der of acceleration is present also in the galactic rotation 
curves which astonishingly successfully can be modeled 
using a modified Newtonian dynamics involving an ac¬ 
celeration parameter umond which again is of the order 
of 10“^ m/s^. Because of this mysterious coincidence of 
characteristic accelerations appearing at different scales 
and due to the fact that all these phenomena appear in 
a weak gravity or weak acceleration regime, it might be 
not clear whether current approximation schemes hold. 
This is one motivation to try to solve the equations of 
motion of test particles in space-times with cosmological 
constant analytically. 

Eurt her more, by looking at the effective potential of 
a point particle moving in the Schwarzschild-de Sitter 
space-time it can be seen that for a certain range of or¬ 
bital parameters a “switching on” of the cosmological 
constant may result in a dramatic change of the orbital 
shape: for a positive cosmological constant bound orbits 
may become escape orbits and for a negative cosmologi¬ 
cal constant escape orbits will become bound orbits, see 
Eig.O The characteristic distance where this happens is 
given by A“^/^ which is of the order of 5 Gpc which is 
roughly the radius of the visible universe and, thus, far 
outside the solar system and our galaxy [8|. However, 
an orbit which is near to the separatrix of Schwarzschild 
geodesics may have a larger sensitivity to a cosmolog¬ 
ical constant which perhaps may not be accounted for 
to the required accuracy in a perturbative approach. In 
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other words, it might be that a comparatively large ac¬ 
celeration cH ^ ^Pioneer ^ <^MOND ^t solar systcm or 
galactic distances may be the result of a very small cos¬ 
mological constant. Therefore, a definite answer to this 
question can be given with the help of an analytical so¬ 
lution only. In addition, the orbits of the Pioneer space¬ 
crafts had been reconstructed using orbit determination 
programs relying on the first order post-Newtonian ap¬ 
proximation. The difference between this approximation 
and the exact orbits with cosmological constant may be 
even more pronounced. 

There is further interest to understand explicitly the 
structure of geodesics in the background of black holes 
in anti-de Sitter space in the context of string theory 
and the AdS/CFT correspondence. In addition, recently 
there also has been a lot of work dealing with geodesics 
and integrability in black hole backgrounds in higher di¬ 
mensions in the presence of a cosmological constant [9|- 

Besides these physically motivated reasons, it is also 
of mathematical importance to derive an explicit analyt¬ 
ical solution of the geodesic equation in a Schwarzschild- 
de Sitter space-time. Orbits of particles and light rays 
have long been used to discuss the properties of solutions 
of Einsteins field equations. In fact, the observation of 
light and particles is the only way to explore the gravi¬ 
tational field. All solutions of the geodesic equation in a 
Schwarzschild gravitational field have been presented in 
a seminal paper of Hagihara M- The solution is given 
in terms of the Weierstrass p-function. With the same 
mathematical tools one can solve the geodesic equation 
in a Reissner-Nordstrom space-time [l^- The analytic 
solutions of the geodesic equation in a Kerr and Kerr- 
Newman space-time have also been given analytically 
(see [m for a survey). Here we expose for the first time 
the complete elaboration of the analytic solution of a 
point particle moving in a Schwarzschild-(anti) de Sit¬ 
ter space-time presented in [T^. Also the entire set of 
possible solutions is described and characterized. For a 
specialized case orbits in a Schwarzschild-(anti) de Sitter 
space-time have been presented Ea. 

Here we consider the general case of geodesics in the 
gravitational field of a spherically symmetric mass in a 
universe with cosmological constant A (of any value), de¬ 
scribed by the Schwarzschild-(anti) de Sitter space-time. 
Because of the static metric and the spherical symme¬ 
try of the problem, the geodesic equation reduces to one 
ordinary differential equation which can be integrated 
formally by means of a hyperelliptic integral. Here we 
explicitly solve this integral. Our calculations are based 
on the mathematically very interesting inversion problem 
of hyperelliptic Abelian integrals studied first by Jacobi, 
Abel, Riemann, Weierstrass, and Baker in the 19th cen¬ 
tury [H-[2l| . The _general ansatz was stated by Kraniotis 
and Whitehouse [22| and Drociuk [ 2 ^ (see also H), but 
in addition to these considerations we explicitly solve the 
equations of motion by restricting the problem to the set 
of zeros of the theta function, the so-called theta divi¬ 



FIG. 1. The tt-component of the Schwarzschild-de Sitter 
metric for various values for A. The dotted line corresponds 
to the Schwarzschild metric. For 0 < A < 1/(9 ?ti^) there are 
two horizons. The dashed line corresponds to the extremal 
Schwarzschild-de Sitter space-time where the two horizons 
coincide. For r < r_ and r > r+ the radial coordinate be¬ 
comes timelike. 


sor. This procedure makes it possible to obtain a one- 
parameter solution of the, in our case, two-parameter in¬ 
version problem. This procedure was suggested by Enol- 
skii, Pronine, and Richter who applied this method 
to the problem of the double pendulum. The resulting 
orbits are classified in terms of the energy and the angu¬ 
lar momentum of the test particle as well as of the value 
of the cosmological constant. A detailed discussion of the 
resulting orbits is given. 

The found analytical solution then is applied to the 
question whether the cosmological constant might be the 
origin of the anomalous acceleration of the Pioneer space¬ 
craft. Over the whole mission, the influence of the cos¬ 
mological constant leads to a modification in the or¬ 
bit of the Pioneers of the order of 10“^m only. The 
found solution is also used to derive the exact post- 
Schwarzschild approximation of the periastron shift. We 
also give one example for the application of this method 
to analytically solve the geodesic equation in higher di¬ 
mensional Schwarzschild, Schwarzschild-(anti-)de Sitter 
or Reissner-Nordstrom-(anti-)de Sitter space-times. 


II. THE GEODESIC EQUATION 


We consider the geodesic equation 

^ cPx^ , „ , dx^ dx^ 


ds‘^ 


( 2 ) 


where ds‘^ = g^^dx^dx^ is the proper time along the 
geodesics and 

{p^} = ^ 9 ^" {dp9au + dagpu - d^gpa) (3) 

is the Christoffel symbol, in a space-time given by the 
metric 


rs 


r 

— r^{d0‘^ + sin^ Od^p), 


ds^ = 1- dt^ — 1-Ar 


rs 


1 


dr^ 

(4) 
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which describes the spherically symmetric vacuum solu¬ 
tion of (HD. This Schwarzschild-de Sitter metric is char¬ 
acterized by the Schwarzschild-radius rs = 2 M related 
to the mass M of the gravitating body, and the cosmo¬ 
logical constant A (unless stated otherwise we use units 
where c = G = 1). The main features of this metric de¬ 
pending on the value of the cosmological constant A are 
shown in Fi g. [J For a general discussion of this metric, 
see e.g. [^,[^. The geodesic equation has to be supple¬ 
mented by the normalization condition = e 

where for massive particles e = 1 and for light e = 0 . 

Because of the spherical symmetry we can restrict our 
consideration to the equatorial plane. Furthermore, due 
to the conserved energy and angular momentum 


^ dt 

E = = 

ds 


( 5 ) 

ds ’ 


( 6 ) 


the geodesic equation reduces to one ordinary differential 
equation 



( 7 ) 

Together with energy and angular momentum conserva¬ 
tion we obtain the corresponding equations for r as func¬ 
tions of s and t 



Equations give a complete description of the dy¬ 

namics. 

Equation ([ 8 ]) suggests the introduction of an effective 
potential 




FIG. 2. The effective potential of a point particle with some 
given L in a Schwarzschild-de Sitter space-time for different 
cosmological constants. 


We rewrite m as 



= P5{U) 


(13) 


with 


P^i^) := —u^ -\- eXu^ + (A(/i — e) p) eXp . (14) 


If not stated otherwise, we take e = 1 in the following. 
Note that p > 0 and A > 0. 

A separation of variables in m yields 


(f- (fO 



u' du' 


(15) 


where uq = u{(Pq). In solving integral ([15]) there are two 
major issues which have to be addressed. First, the in¬ 
tegrand is not well defined in the complex plane because 
of the two branches of the square root. Second, the solu¬ 
tion u{ip) should not depend on the integration path. If 
7 denotes some closed integration path and 


C udu 

^ \/ Pbi^P) 


(16) 


this means that 


shown in Fig. [21 It is worthwhile to note that for light, 
i.e. e = 0 , the cosmological constant just gives a constant 
contribution to the effective potential and, thus, does not 
influence m and 0 . However, it still influences the 
motion of light through the timing formula ([9]) . 

As usual, we introduce a new variable u = rs/r and 
obtain 

A ^ 1 

= eXu^ {X{p - e) + p) ^ eXp^ (11) 

J 

with the dimensionless parameters 

A ^ , n:= and p := ^Arl. ( 12 ) 



ip-ipo-u> = 



u' du' 

Twy 


(17) 


should be valid, too. Hence, the solution u{(p) of our 
problem has to fulfill 


u{p) = u{(p — uj) (18) 

for every cj 7 ^ 0 obtained from an integration (HSl). A 
function u with the property ([T 8 |) is called a periodic 
function with period uj. These two issues can be solved 
if we consider Eq. m to be defined on the Riemann 
surface X of the algebraic function x yj P 5 (x). 
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III. THE INVERSION PROBLEM 


Let X be the compac t Riem annian surface of the al¬ 
gebraic function x ^ P 5 (x). It can be represented as 
the algebraic curve 


X:={z = (x, y) e C2 12/2 = ( 19 ) 


( 2 ^ or as the analytic continuation of VA- The last one 
can be realized as a two-sheeted covering of the Riemann 
sphere which can be constructed in the following way: 
let e^, i = 1,..., 5, be the zeros of P 5 and ee = 00 (for a 
polynomial of 6 ^^ order the zero is finite). These are 
the so-called branch points. Now take two copies of the 
Riemann sphere, one for each of the two possible values 
of \/^, and cut them between every two of the branch 
points ei in such a way that the cuts do not touch each 
other. These are the so-called branch cuts, see Fig.[3l Of 
course, the two copies have to be identified at the branch 
points where the two values of identical. They 

are then glued together along the branch cuts in such a 
way that \/^ together with all its analytic continuations 
is uniq uely de fined on the whole surface. On this surface 
X y^P^{x) is now a single-valued function. This con¬ 
struction can be visualized as a ’’pretzel”, see Fig.[3l For 
a strict mathematical description of the construction of 
a compact Riemannian surface, see [29|, for example. 

Every compact Riemannian surface can be equipped 
with a homology basis {a^, bi\i = 1,..., ^} G Hi{X, Z) 
of closed paths as shown in Fig. [3l where g is the genus 
of the Riemannian surface, see the next section. From 
the construction of the Riemannian surface it is already 
clear that integrals over these closed paths indeed do not 
evaluate to zero and, hence, have to be periods of the 
solution of m- The task now is to analyze the details 
of periodic functions on such Riemannian surfaces. 


A. Preliminaries 


Compact Riemannian surfaces are characterized by 
their genus g. This can be defined as the dimension of 
the space of holomorphic differentials on the Riemannian 
surface or, topologically seen, as the number of ‘holes’ of 
the Riemannian surface. Let Pd = ^ Poly¬ 

nomial of degree d with only simple zeros and X be the 
Riemann surface of y/Pd. Then the genus of X is equal 
to ^ , where [x] denotes the largest integer less 

or equal than x 0 . Hence, in our case of P 5 the genus 
of the Riemann surface is g = 2. 

In order to construct periodic functions on a Riemann 
surface we first have to define a canonical basis of the 
space of holomorphic differentials {dzi\i = 1 ,...,^} 
and of associated meromorphic differentials {dri\i = 




FIG. 3. Riemannian surface of genus g — 2^ with real branch 
points ei,...,e6. Upper figure: Two copies of the complex 
plane with closed paths giving a homology basis {ai^hi \ i — 

1.. .. ,p}. The branch cuts (thick solid lines) are chosen from 
e 2 i-i to e 2 z, i = 1,...,^ + 1. Lower figure: The ’’pretzel” 
with the topologically equivalent homology basis. 

1.. .., ^} on the Riemann surface by 

^ x'^~^dx 

dzi := — , , 

with Xj being the coefficients of the polynomial Pd 
In our case these differentials are given by 

_ dx _ xdx 

dzi := , , dz 2 := , , 

^ 3x^ — 2x‘^ + Ax _ _ x‘^dx 

dvi '.= - . - dx , dT 2 := - , , 

where A is defined in ([12]). We also introduce the period 
matrices (2 uj^2uj') and (2g^2g') related to the homology 
basis 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 




(24) 


The differentials in (|2Q|) and m have been chosen such 
that the components of their period matrices fulfill the 
Legendre relation 


fuj uj'\ f 0 —tq\ f 00 cj'V 1 . f 0 — 

[g Vj [h 0 ) [g g') “ “ 2 ^' [ig 0 j ’ 

(25) 
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where Ig is the g x g unit matrix, j^ . 

Finally we introduce the normalized holomorphic dif¬ 
ferentials 


dv := {2uj) 


( dzi\ 

dZ2 


\dzgj 


(26) 


The period matrix of these differentials is given by 
where r is defined by 


r \= uj . (27) 

It can be shown [3l| that this normalized matrix r always 
is a Riemannian matrix, that is, r is symmetric and its 
imaginary part Imr is positive definite. 


B. Jacobi’s inversion problem 


Let us consider now the Abel map 

Axq : X Jac(X), X f dz (28) 

J Xo 

from the Riemannian surface X to the Jacobian 
Jac(X) = C^/r of X, where F = {ujv-^uj'v' | G 
is the lattice of periods of the differential dz. The image 
Axo (X) of X by this continuous function is of complex di¬ 
mension one and, thus, an inverse map A~^ is not defined 
for all points of Jac(X). However, the ^-dimensional 
Abel map 


9 pxi 

: S^X Jac(X), (rcj,..., Xgf 

i=l 


(29) 

from the ^fth symmetric product S^X of X (the set of 
unordered ’’vectors” (xi, ... ^XgY where Xi G X) to the 
Jacobian is one-to-one almost everywhere. Jacobi’s in¬ 
version problem is now to determine x for given (f from 
the equation 


0=A^^{x). 


(30) 


In our case g = 2 this reads 
dz 




^2 = 


-I. 


VP^iz] 

zdz 




f 

Jxo 

rx2 

I 

9 Xn 


dz 


zdz 


(31) 


We will see later, that we can solve our problem ([T5]) as 
a limiting case of this Jacobi inversion problem. 


C. Theta Functions 

The Riemannian surface of genus g has 2g independent 
closed paths, each corresponding to a period of the func¬ 
tions defined on these surfaces and, hence, to a period of 


the solution u of ([T5]) . In order to construct 2^-periodic 
functions, we need the theta function'd : ^ C, 


'd{z;T) := 

rhe'La 


{Trh-\-2z) 


(32) 


The series on the right-hand side converges absolutely 
and uniformly on compact sets in (D^ and, thus, defines a 
holomorphic function in (D^. This is obvious from the es¬ 
timate Re(m^(iT)m^) < —cm^m for some constant c > 0, 
what follows from the fact that Re(ir) is negative defi¬ 
nite. The theta function is already periodic with respect 
to the columns of \g and quasiperiodic with respect to 
the columns of r, i.e., for any n e the relations 

'd{z-\- IgH; r) = d{z] r) , (33) 

diz + Tfi-, t) = r) (34) 


hold. We will also need the theta function with 
characteristics [ 32 I g^h E defined by 


d[g,h]{z;T) := Y, 


rhE'Ld 




Later it will be important that for every g, h the set 
G O I 'd[g^h]{z]T) = 0}, Called a theta 
divisor^ is a {g — l)-dimensional subset of Jac(X), see 
[m or (j49j). 

The solution of Jacobi’s inversion problem (|3T]) can be 
explicitly formulated in terms of functions closely related 
to the theta function. First, consider the Riemann theta 
function 


de{x;r):=dyj dv-e^rj, (36) 

with some arbitrary but fixed e G C^. The Riemann 
vanishing theorem, see e.g. [3l|, states that the Riemann 
theta function is either identically to zero or has exactly 
g zeros xi,..., for which 

Y! dv = e + K^^ (37) 

i=l '^^0 

holds (modulo periods). Here Kxq G (D^ is the vector 
of Riemann constants with respect to the base point xq 
given by {rjj is the jth diagonal element) 

Kxo,j = “ Z / ([ 

^ l^j ^dxo / 

If the base point xq is equal to oo, this vector can be 
determined by 


Koc 



(39) 
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where e 2 i is the starting point of one of the branch cuts 
not containing oc for each see [s^. Hence, K^o can be 
expressed as a linear combination of half-periods in this 
case. For problems of hyperelliptic nature it is usually 
assumed that the Riemann theta function 'de does not 
vanish identically. However, here we are interested in 
the opposite case: we want to restrict Jacobi’s inversion 
problem (|3T|) to the set of zeros of 'd(' -h Kxq]t)^ which is 
called the theta divisor 0^ . 

The solution of ()3ip and, thus, of m can be formu¬ 
lated in terms of the derivatives of the Kleinian sigma 
function cr : ^ C, 

a{z) = ((2c^)-i£+ , (40) 

where the constant C can be given explicitly, see jsoj . 
but does not matter here. Jacobi’s inversion problem can 
be solved in terms of the second logarithmic derivative 
of the sigma function called the generalized Weierstrass 
functions 


d d 


(Ji{z)crj{z) - a{z)aij{z) 




0-2 (i) 


(41) 

where denotes the derivative of the sigma function 
with respect to the i—th component. 


D. The solution of the Jacobi inversion problem 

The solution of Jacobi’s inversion problem (|3Q|) can be 
given in terms of generalized Weierstrass functions. Let 
X be the Riemannian surface of ^/P where P is defined 
by P{x) := (this can always be achieved by 

a rational transformation). Then the components of the 
solution vector x = (xi,..., XgY are given by the g solu¬ 
tions of 

-^2:5 - = 0 , (42) 

i=l 

where (f is the left hand side of (|3Q|) . Since x e S‘^X 
there is no way to define an order of the components of 
X. In our case of = 2, we can rewrite this result with 
the help of the theorems by Vieta in the form 

Xl PX2 =^p22((p) , 

' 4 («) 

XiX2=- —pl2{(f) . 

^5 

IV. SOLUTION OF THE EQUATION OF 
MOTION IN SCHWARZSCHILD-(ANTI-)DE 
SITTER SPACE-TIME 

Now we apply the results of the preceding section to 
the problem of the equation of motion in Schwarzschild- 
(anti-)de Sitter space-time, Eqs. ([T3|) and ([T5]) . As al¬ 
ready mentioned before, the solution of the equation of 


motion can be found as a limiting case of the solution of 
Jacobi’s inversion problem in the case of genus g = 2. 


A. The analytic expression 


To start with, we rewrite Jacobi’s inversion problem 
(|3T]) in the form 


02 


L 

ioo \J 


’ oo y Pbix) 
xdx 


+ 


+ 


r" dx 

L 

xdx 

loo ’ 


(44) 


where 


P- 


2 [ dz. 
Juq 


(45) 


Note that the right-hand side of (|44|) is exactly Aoo{u), 
the image of the Abel map defined in (|29l). i.e. 0 = 
Aoo{u). We use the obvious identity (compare [25[) 


Ui 


y UlU2 

hm - 

U2^oo Ui + U2 


(46) 


and insert the solution of Jacobi’s inversion problem (liSll . 
Then 


Ui 


hm 

P22(<(>) 


= lim 

U2^00 


a 


(0)0-12 (y - O-i(0)o-2(0) 


- 0-0-22(0) 


where 


’-(0oo)o'l2(0 OO ) - O'l(0 OO )O'2(0 o 
O-i(0oo) - O-(0 OO )cr22(0 oo) 


= lim (j) = Aoo(uoo) 

U2^00 


(47) 


(48) 


with Uoo = (^). Note that the definition of the sigma 
function (HUl) and, hence, of the generalized Weierstrass 
functions dm includes the vector of Riemann constant 
Kxq with xo = oo in our case, which is given by Koo = 

^{ 1 / 2 ) + (1/2) (see dm or (SH). 

The above limiting process also transfers Jacobi’s 


inversion problem to the theta divisor 0 


Ko. 


With 


{2uj) = {2uj) ^Aoo{uoo) = dv and the theorem 


id 




_(v0 41/2)] {z-,t)=0 

r 

3x:z= dv 

J 00 

proven by Mumford [3l| it follows that 

0 = ^ (1/0’(1/2) ((2yVoo;'r). 


(49) 


(50) 
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Via Eq. (j35j) this is equivalent to 


0 = ^ 


+ (i? 2 ) 


((2w) 

and with (gni) this means 

cr(0oo) = 0. 

We first use this result in (|T7|) and obtain 

( 0 oo) 


Ui 


(^2{4>oo) 


(51) 


(52) 


(53) 


Theorem (Hi) also tells us that (2uj) 

i.e 


-1 


(j)oo is an ele¬ 


ment of the theta divisor 0^^, i-e. the set of zeros of 

'd (^ 1 / 2 ) 5 ( 1 / 2 ) : that, in the case ^ = 2, 0^^ is 

a manifold of complex dimension one. Note that the re¬ 
striction to the theta divisor is only possible because 00 
is a branch point what is essential for the validity of the¬ 
orem (|49|) . 

Since 0^ is a one-dimensional subset of (D^, there 
is a one-to-one functional relation between the first and 
the second component of (2co’)“^0oo- By the definition of 
0OO in (07]) and Eq. (05|) we have 

0OO = lim 0 

U2^00 


r 

= lim 0—21 

U2^00 

J Un J U 


dz 


Uo 

00 


dz. 


'uo Juq 

The physical coordinate (f is given by (USD, 
zdz 




-f 

J Uo 


-L 


^ Uo \/Pb (^) 

We insert this in (|54|) and obtain 


+ ^0 — / + ^0 • 


(54) 


(55) 


dzi — dzi 

Juo ^ ^ 

^ - ^0 - dZ2 


Juo 

rOO 


Cdzi-J-dz^ 
Po 


P- 


(56) 


where (p^ = dz 2 depends only on the initial values 

uq and cpo. We choose for each (p a cpi such that 0q := 

( ) is equal to 0oo- Then {2uj)~^0q = {2uj)~^$oo 

yP ~ 

is an element of the theta divisor 0^ and we finally 
obtain 


r{p) = 


rs O'2((?oo) 0-2 (i^e) 

- = —rg- ^ - = “"Ts- • 

'^ip) cfi{4>oc) cri((^e) 


(57) 


This is the analytic solution of the equation of motion 
of a point particle in a Schwarzschild-(anti-)de Sitter 
space-time. This solution is valid in all regions of the 
Schwarzschild-(anti-)de Sitter space-time and for both 
signs of the cosmological constant and can be computed 
with arbitrary accuracy. The explicit computation of the 
solution is described in Appendix A. 


B. Light trajectories 


In the case of light trajectories, the situation simplifies 
considerably. The equation of motion is then given by 

(^) = + V + =;-P3(w) • (58) 

Light rays are uniquely given and, thus, uniquely charac¬ 
terized by the extremal distance to the gravitating body, 
that is, the smallest or largest distance (in the case of 
a Schwarzschild space-time, it is also possible - due to 
its asymptotic flatness - to take the impact parameter as 
characteristic of a light ray). This extremal distance uq 
is characterized by 


which gives 
of motion is 


du 

dp 


= 0 , 


U=Uo 


(59) 


+ A/i + = 0. Then our equation 


=u^ -V? -ul^ul (60) 

dp J 


which is the same type of equation as in Schwarzschild 
geometry. With a substitution = 4x -h | this reads 



are the Weierstrass invariants. This differential equation 
can be solved directly in terms of elliptic functions, i.e. 


r{p) = 


rs 




rs 


u{p) Ax{p) -h \ 


4p(^-^0;^2,^3) + ^ ’ 
(63) 

where p is the Weierstrass function [33|, and p'^ 
is given by the initial values po and xq, Po = Po + 
The corresponding light trajectories have 

been exhaustively discussed in [H. 

In a recent paper [35| Rindler and Ishak discussed 
light deflection in a Schwarzschild-de Sitter space-time. 
Though the equation of motion is the same as in 
Schwarzschild space-time, the measuring process for an¬ 
gles reintroduces the cosmological constant in the ob¬ 
servables. According to their scheme, the exact angle 
between the radial direction and the spatial direction of 
the light ray is now given by 


tan 7 /; = 


1 - #5 - 


\l^(> 


IS ^ _ fl _ r^\ 

ro J \ r{(p) J 


(64) 


where in the expression dr j dp from ([7|) the ^ -h ^A has 
been replaced by the ro related to uq . This now is valid 
for all light rays, not only for those rays showing a small 
deflection as discussed in j35| . 
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V. THE CLASSIFICATION OF THE 
SOLUTIONS 

A. General classification 

Having an analytical solution at hand we can explore 
the set of all possible solutions in a systematic manner. 
The shape of an orbit depends on the energy E and the 
angular momentum L of the particle under consideration 
as well as the cosmological constant A (the Schwarzschild 
radius has been absorbed through a rescaling of the ra¬ 
dial coordinate). These quantities are all contained in 
the polynomial Ps{u) through the parameters A, /r and 
p m- Since r (and u) should be real and positive it is 
clear that the physically acceptable regions are given by 
those u for which E > Ve^. The zeros of P 5 are related 
to the points of intersection of E and , and a real and 
positive P 5 is equivalent to P > as can also be seen 
from ([13]). Hence, the number of positive real zeros of P5 
uniquely characterizes the form of the resulting orbit. 

Since P5 goes to —00 if x ^ —00 and to oc if x ^ 00, 
P5(0) is positive if the number of positive real zeros of P5 
is even and negative if it is odd. If we denote by ei,..., 
the positive real zeros, then it follows that the physically 
acceptable regions are given by [ 0 , ei], [02, 63], ..., [e^, 00] 
if n is even and by [ei, 62],... , [e^, 00] if n is odd. With 
respect to r we have the following classes of orbits (see 
Fig. gj): 

(i) the region [ 0 ,ei] corresponds to escape orbits, 

(ii) the region [e^,oo] corresponds to orbits falling into 
the singularity, i.e. to terminating orbits, and 

(hi) the regions [ei,e^+i] correspond to bound orbits. 

This means that for any arrangement of zeros of P5 there 
exist terminating orbits. Furthermore, for an even num¬ 
ber of positive real zeros we have escape orbits and for 
more than three real positive zeros we have bound orbits. 
The case that there is no positive real zero corresponds 
to a particle coming from infinity and falling to the sin¬ 
gularity, see Fig. gJ Quasiperiodic bound orbits exists 
only if there are three or more positive zeros. 

It can be shown that there are no more than four real 
positive zeros for our polynomial (HI: We decompose the 
polynomial P5 into its (in general complex) zeros P5 (u) = 

{u — ui){u — U 2 ){u — us){u — U 4 ){u — Us). Multiplication 
and comparison of the coefficients of the terms linear in 
u yields 

UiU2UsU4EUiU2UsUsEUiU2U4Us-\-UiUsU4UsEU2UsU4Us = 0 . 

(65) 

The assumption that all zeros are real and positive con¬ 
tradicts Eq. ([65]) . Therefore, in any case there are at 
most four real positive zeros. 

Figure [5] shows the arrangement of zeros of Ps{u) for 
some chosen values of p as a (p. A) - diagram. The code of 
grayscales is as follows: black corresponds to four, dark 
gray to three, gray to two, light gray to one and white to 
no positive real zero. 


P 5 - 


P5‘ 



u 




(a) No real positive zero (b) One real positive zero 



(c) Two real positive zeros (d) Three real positive 

zeros 



(e) Four real positive 
zeros 


FIG. 4. The five possibilities of having real positive zeros of 
Psiu). The allowed regions of particle motion are shaded in 
gray. The zeros correspond to the zeros of Uff = E in Fig. [2] 
(note that = 0 corresponds to r = oo and = oc to r = 0). 
Bound nonterminating, quasiperiodic orbits exist only if there 
are three or more positive zeros. 


B. Discussion with respect to A 

Based on Fig.[5]we are now discussing the orbits related 
to different values of A. Each Plot in Eig. [5] comprises 
all effective potentials (of the form shown in Eig. [2|) for 
all possible values of L and all particle energies E and, 
thus, contains the complete information about all orbits 
in Schwarzschild-(anti-)de Sitter space-times for a given 
value of A. 

a) Let us first consider A = 0. In this case the con¬ 
stant term in Ps{u) vanishes and Ps{u) = u‘^Ps{u) 
so that rt = 0 is a zero of Ps{u) with multiplicity 2 . 
The polynomial P 3 corresponds to the Schwarzschild 
space-time and has been extensively discussed in H- 
Nevertheless, let us examine some regions of (A, ja) 
and possible orbits so that they can be directly com¬ 
pared with orbits for A 7 ^ 0 . As seen in Eig. EKd), the 
straight line /i = 1 divides the plot in two parts. Eor 
/i < 1 there is an odd number of positive real zeros, 
i.e. we may not have any escape orbits in this regions. 
The light gray area corresponds to one real positive 
zero and, therefore, it is only a terminating orbit pos¬ 
sible whereas in the dark gray region there may be in 
addition a bound orbit. Eor /i > 1 there is an even 
number of positive real zeros and, thus, there is al- 
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FIG. 5. The zeros of P^iu) in a (/x, A)-diagram for different values for the cosmological constant (/x is along the x-axis, A along 
the y-axis). The gray scales encode the numbers of positive real zeros of the polynomial P 5 : black = 4, dark gray = 3, gray = 
2 , light gray = 1, white = 0 . In the plot for A = 0 characteristic lines are shown (the left upper point of the dark gray region 
has the coordinates (/x = f ? ^ = f); fhe upper intersection point of the /x = 1 line with the dark gray region is at A = 


ways an orbit which reaches infinity. The gray region 
corresponds to two real positive zeros and, hence, to 
a terminating and an escape orbit. The white region 
represents the case of no positive real zeros, i.e. an 
orbit which comes from infinity and falls into the sin¬ 
gularity. Also, beside xx = 0 there is a further real zero 


with multiplicity 2 on the straight line /x = 1. 

b) Let us compare now this with the case A > 0, see 
Fig. [5fe)-(g). We immediately recognize that left to 
the /X = 1 line the plot significantly changed. In ad¬ 
dition, we notice that for growing A this straight line 
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shifts a bit to the left. Left of the straight /i-line there 
is now one more positive real zero in each region. This 
means that a particle which for A = 0 is in a light gray 
region now is in a gray region and, thus, may reach 
infinity. The same happens in the region which was 
dark gray for A = 0 and is now black. A particle with 
/i < 1 in the small band now right of our straight line 
switched to a gray or white region depending on its A 
value. 

For a large positive cosmological constant the black 
area will disappear, that is, there will no longer be any 
bound orbit. This is clear from the following: First we 
introduce f := r/rs- For A = 0 the effective potential 
lAeff possesses two different extrema f± if A < 1 The 
smaller extrema f_ is a maximum whereas the larger 
is a minimum. The extrema are bounded by f_ <3 
and f+ > I — 3 > 3. If we add the term containing 
A > 0, which is of the form of a parabola, f_ shifts to 
the left and to the right. Thus, f_ < 3 and > 3 
remain valid. In general, a second maximum ta > 
will appear and — oc for r ^ oo. It follows that 

there will be no bound orbit if the minimum and, 
thus, the maximum va disappears. This is fulfilled if 
we choose A such large that the gradient is negative 
for all f > 3, for example Ar| > ~ Since 

^ > 3, the choice Ar| > ^ ensures that for any choice 
of A < -^ there will be no bound orbits. This is of 
course only a rough estimate which can be improved. 

c) If A < 0 , the situation changed the other way around. 
We again immediately see that the right side of the 
plot significantly changed. The straight /i-line is no 
longer so easy to identify, but if we take into account 
the number of all real zeros, we can say that it shifts to 
the right when the absolute value of A growths. An 
exception to this is the part for small A. There the 
line bends to the right and allows a switch from the 
gray part for A = 0 to the light gray part of A < 0 . 
Nevertheless, we can say that on the right side of the 
plot we have now an additional real positive zero and, 
thus, also an odd number of positive real zeros. This 
means that a particle is no longer able to reach infinity 
for any (A,/i). In the region being white for A = 0 
and which now is light gray we have now a bound 
terminating orbit. In the for A = 0 gray region which 
now is dark gray the escape orbit becomes bound. 


C. Plots of orbits 

In Figs. [6]l8] some of the possible orbits are plotted. The 
figures are organized in order to highlight the infiuence 
of the cosmological constant. For all orbits in each figure 
the parameters /r and A, that is, E and L, are the same. 
The absolute value of the cosmological constant is chosen 
as |A| = 10“^ in all plots. All plots are created from the 
analytical solution derived in Sec. IIV A1 


In Fig. [ 6 ] the parameters are /i = 0.92 and A = 0.28 
which belong to the dark gray region in Fig. [5fd). For 
a vanishing cosmological constant this corresponds to a 
bound periodic orbit and to a bound terminating or¬ 
bit ending in the singularity. The corresponding orbits 
are shown in Fig. [ 6 fa-b). For a positive cosmological 
constant the overall structure changes considerably since 
there will be a third type of orbits not present in the 
Schwarzschild case. Beside the terminating and bound 
orbits in Fig. me) which both look quite similar to the 
corresponding orbits in the Schwarzschild case there is an 
escape orbit which is repelled from the potential barrier 
related to the positive cosmological constant. Fig. [D^e). 

The next parameter choice is fi = 1.1 and A = 0 . 2 . 
For vanishing A these parameters lay in the gray area 
of Fig. [SJd) denoting two zeros and, thus, correspond to 
a quasihyperbolic escape orbit. Fig. [7fa), and a termi¬ 
nating orbit ending in the singularity. Fig. [7fb). For the 
chosen A = —10“^ km“^ the situation changes dramati¬ 
cally as can be read off from Fig. [5fa) in comparison to 
Fig.E^d): Now we have three zeros and, thus, one bound 
orbit and one terminating orbit ending in the singularity, 
see Figs. [TJc-d). Switching on the negative cosmologi¬ 
cal constant makes an escape orbit a bound orbit. This 
of course has to be expected as one can see from Fig. [2] 
that there are no escape orbits for negative cosmological 
constant. 

Our third choice of parameters is /r = 0.8 and A = 0 . 2 . 
For A = 0, Fig. [5fd) this lays in the light gray region of 
one zero where is a terminating orbit only. This orbit is 
shown in Fig. [ 8 fa). For a positive cosmological constant 
A > 0 these parameters are in a gray region with two 
zeros indicating a terminating and an escape orbit, see 
Figs. [Hl^b-c). The orbit in Fig. [ 8 l^c) again is a reffection 
at the A-barrier. 


VI. ON THE PIONEER ANOMALY 

We apply the obtained analytical solution in order 
to decide whether a nonvanishing cosmological constant 
may have an observable inffuence on the Pioneer satel¬ 
lites. From we may deduce the energy and an¬ 
gular momentum of the Pioneers after their last ffybys 
at Jupiter and Saturn, respectively, with respect to the 
barycenter of the inner solar system, i.e. the Sun, Mer¬ 
cury, Venus and Earth-Moon. This means that we used 
the value 

2 GM , , , 

rs = —^ = 2.953 266 762 363 45 km, 66 ) 

for the Schwarzschild radius, derived from GM = 

I . 000 005 65/c^(AU^/day^) with Gauss’ constant k = 
0.017 202 098 95 defining the astronomical unit AU. Here 
all numbers are taken with 12 digits what corresponds to 
the today’s accuracy of solar system ephemerides. 

In the case of Pioneer 10, the velocity at infinity = 

II. 322 km s~^ taken from gives us the energy per 
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(a) A = 0 , ro = 20.951km (b) A = 0 , ro = 5.010km 



FIG. 6. Orbits for fi = 0.92 and A = 0.28. The upper row is for vanishing, the lower row for positive A. (a) and (c): bound 
orbits with perihelion shift, (b) and (d) terminating orbit ending in the singularity, (e): Reflection at the A-barrier. There is 
no analogue of (e) for A = 0. Black circles always indicate the Schwarzschild radius. 


unit mass Em = H—and therefore the parameter 

/i, 2 

E^ 

= ^ = 1.000 000 00143. (67) 

The angular momentum per unit mass is given by Lm = 

gu, where v = ^2GM Q + is the velocity at peri- 

apsis distance g; a is the semimajor axis. From this we 
derive the parameter A, 

^2 2 

A = ^ = 2. 855 572 373 82 • 10“^ . (68) 

In the case of Pioneer 11 we obtain for the parameters fi 
and A 

li = 1.000 000 00122, 

( 69 ) 

A = 1.340 740 574 59-10“^ ^ ^ 

With these coefficients we now can determine the ex¬ 
act orbits of Pioneer 10 and 11 in the cases A = 0 and 


A = 10“^^km“^. From these exact orbits we calculated 
the differences in position (in m) for a given angle (f (in 
rad) and the difference in the angle (in rad) for a given 
distance r (in m) of a test particle moving in a space- 
time with and without cosmological constant. The Pi¬ 
oneer anomaly appeared in a heliocentric distance from 
about 20 to 70 AU. For r in this range, we compute now 
the difference (Pa=o{t) — ipAy^o{r) in azimuthal position 
with and without cosmological constant for both craft. 
Regarding Pioneer 10, the difference is in the scale of 
10“^^rad, which corresponds to an azimuthal difference 
in position of about 10“^m. For Pioneer 11, the differ¬ 
ence is in the scale 10“^^rad, which corresponds to an 
azimuthal difference in position of about 10“^m. 

The range of 20 to 70 AU corresponds to an angle 
between 0.47r and O.Gtt if ^po = 0 corresponds to the pe- 
riapsis. In this range, we compute the radial difference 
rA=o{^) —'rA^o(^) also for both craft. For Pioneer 10 we 
obtain a difference in the scale of 10“^m, for Pioneer 11 
in the scale of 10“^m. 

Therefore we can say, that for the present value of the 
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(a) A = 0, ro = 6.776km 


(b) A = 0, ro = 4.642km 



(c) A = -10-5km-2^ ro = 185.37km (d) A = -IQ-^km-^^ ^*0 = 4.639km 

FIG. 7. Orbits for /x = 1.1 and A = 0.2. The upper row is for vanishing, the lower row for negative A. 


cosmological constant the form of the Pioneer 10 orbit 
practically does not change. For a definite estimate of 
the differences of the Pioneer orbits in Schwarzschild and 
Schwarzschild-de Sitter space-time one of course has to 
analyze the time course of these orbits. However, the 
time variable is influenced by the cosmological constant 
in the same way as the radial coordinate so that no 
change in our statement will occur. Therefore, the in¬ 
fluence of the cosmological constant on the orbits cannot 
be held responsible for the observed anomalous accelera¬ 
tion of the Pioneer spacecraft. 


VII. PERIASTRON ADVANCE OF BOUND 
ORBITS 

In the case that P 5 has at least three real and positive 
zeros, we may have a bound orbit for some initial values. 
The periastron advance Aperi for such a bound orbit is 
given by the difference of the 27r-periodicity of the angle 
(f and the periodicity of the solution r{(p) (which is the 
same as the periodicity of u{(p)). Let us assume that the 
bound orbit corresponds to the interval [e/c, e/c+i], where 


Ck and e/e+i are real and positive zeros of P 5 , and that the 
path ai surrounds this real interval. Then the periastron 
advance is given by 

r^k+i xdx 

Aperi = 27r - 2 uj2i = 2tt - 2 / , (70) 

Jek VP5{x) 

where 2 uj 2 i is an element of the (canonically chosen) 2x4 
matrix of periods (2co’, 2uj') of \/^, see Eq. (|24|) . We now 
calculate the post-Schwarzschild limit of this periastron 
advance in the case that the considered bound orbit is 
also bound in Schwarzschild space-time. 

For doing so we first expand x / ^/^'(^ to first order 
in A 

^ ^ K 2 + ^ 

yAF) yFF) 6 ®a;2P3(a:)yFF) ’ 

where P 3 (x) = x^ — x^ -j- Ax -h A(ju — 1) is the polynomial 
for the corresponding Schwarzschild case given by A = 0 . 

In the next step we have to integrate both terms involv- 
ing P 3 within the Weierstrass formalism, see for example 
. Employing the substitution x = 4z^l/3we rewrite 
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(a) A = 0, ro = 3.625km (b) A = 10 ^km ro = 3.625km (c) A = 10 ^km ro = 237.61km 

FIG. 8 . Orbits for fi = 0.8 and A = 0 . 2 . The upper graph is for vanishing, the lower graphs for positive A. There is no analogue 
of (c) for A = 0 . 


P3 in a Weierstrass form 


P^ix) = 4?{Az^ - g 2 Z - gs) = i^Pw{z ), (72) 


where 


-l_lx 

12 4 ^ 

“ 16 3"^ ~ 


(73) 

(74) 


are the Weierstrass invariants. We assume that the 
orbit under consideration is bound not only in the 
Schwarzschild-de Sitter but also in the corresponding 
Schwarzschild space-time. This means that the three 
largest real zeros of P5 are positive and, thus, the zeros 
zi > Z 2 > zs > —^ of Pw are all real. The square root 
y/Pw is branched over zi, 2)2 and 2)3 and, thus, the ellip¬ 
tic function p based on y/Pw has a purely real period uji 
and a purely imaginary period uj 2 ‘ They are given by 


CJi 


002 


dz 


A yJPw{z) 


i 


dz 

y^Pw{z) 


(75) 


where the path A runs around the branch cut from 2)3 to 
Z 2 and the path B around Z 2 and 2)1, both clockwise. 
The branch of the square root in (USD is chosen such 
that yJPw > 0 on [2)3,2)2] and, thus, y/Pw negatively 
imaginary on [2)2, Zi]. The branch points of p can be ex¬ 
pressed in terms of the periods: Zi = p(pi), Z 2 = p{p 2 ) 
and 2)3 = p{p 3 ) with pi = uJi/ 2 , p 2 = {uji Puj 2)/2 and 
Ps = 0 J 2 I 2 . The fundamental rectangle in the complex 
plane spanned by the periods uji^uj 2 of p is denoted by 
R = {xuji + yuj 2 I 0 < X, ^ < 1}, see Fig. [9l 

Let the three biggest real and positive zeros of P5 be 
given by > X 2 > X 3 >0. Then, for the canonical 



FIG. 9. The fundamental rectangle R 


choice of the matrix of periods uj of y/P^^ the integration 
path ai runs from X 3 to X 2 and back with conversed sign 
of the square root. Let the path 7 be the preimage of ai 
by 14 i-A p{u) = ^ in the fundamental rectangle R. For a 
positive cosmological constant A, we have X 3 < Z 3 < Z 2 < 
X 2 and, thus, 7 starts at some purely imaginary 7 ( 0 ) = 
i4i G P with 0 < Im(i4i) < Im(p3) and goes straight to 
7 ( 1 ) = ui P uji. Then, for any rational function P, we 
obtain 

/ / F{p{A)du. (76) 

Jai ^yPw[Z) 

This is derived from the differential equation 

= yF7(p(7)), (77) 
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where the branch of the square root was chosen to be 
consistent with the sign of p'. 

The integration of the first part on the right-hand side 
of ([TT]) is straightforward and yields the Schwarzschild 
period 


= 1.1849627128268641 x lO-^^niVs^ for the an- 

gular momentum per unit mass Lm and \/Em = 
0.029979245417779875 x 10^^ m/s for the energy per unit 
mass Em given in [22| . These values lead to the zeros 


/ -^=1 /*■=„,. (78) 

Jai JPs (x) Jai JPw{z) Jui 


The integration of the second part on the right-hand side 
of ([7T]) is more involved and is performed in Appendix 
B. As a result we obtain the first order approximation of 
the periastron shift with respect to A: 


^peri — 27r 


xdx 




= 2^- ^,+A^ 


+ A 


rji + ZjCOi 


2rji - pi 6 p"{uo) 

16p'(mo)^ 16 p'(wo)® 


(4^i + 3)" 


(-r^i^o + C(mo)) 


+ 0{K^), 


(79) 


where uq is such that p(mo) = “H' 

The terms in this expression involving pj and uq can 
partly be replaced by terms containing the Weierstrass 
invariants g 2 and gs. From the differential equation (f77|) 
we derive 


p'K) = \/4p(uo)^ -52p(«o) -53 = Y■ 

(80) 

The first derivative of dZZl) yields 2 p" = 12p^ — ^2 and, 
thus, gives 

p''{pj) =Qz^-^92 and p"(uo) = ^ - ^52 , ( 81 ) 

where the g 2 ^ ^3, as well as the zeros of Pw can be ex¬ 
pressed by p and A. 

The result ([79]) gives the post-Schwarzschild perias¬ 
tron shift in a closed algebraic form. The advantage of 
this result is that no further integration is needed. An¬ 
other advantage lies in the fact that only elliptic func¬ 
tions and related quantities are used which are well de¬ 
scribed and tabulated in mathematical books and which 
are also well implemented in common commercial math 
programs. What is still left to do is to express the 
result ([79]) in terms of, e.g., rmin and rmax or, equiv¬ 
alently, in terms of the semimajor axis and the eccen¬ 
tricity. These quantities are directly observable and also 
have the advantage that an expansion in terms of m/rmin 
and m/rmax can be performed giving in addition a post- 
Newtonian expansion. This will be described elsewhere. 

Let us apply these formulas to the perihelion ad¬ 
vance of Mercury and compare with the results of 
Kraniotis and Whitehouse, We take the val¬ 

ues rs = 2953.25008m for the Schwarzschild-radius, 


= 0.166666640041880, 

Z2 = -0.083333317283501, 

^3 = -0.083333322758379 

of Pw and to the periods 

cji = 3.1415929045225246, 
a;2 = 20.40939163938517992, 

T = 6.49651061090841872, 

which all compare well to the results in [22|. Also the 
physical data, i.e. the aphel the perihel rp and the 
perihelion advance in Schwarzschild-space-time A^, com¬ 
pare well to [22I and also to observations [lH : 

rA = 6.981708938652731 • lO^^m, 
rp = 4.600126052898539 • 10^ V, (82) 

= 42.980165 arcsec cy“^. 

Here we used the rotation period 87.97 days of Mer¬ 
cury and 100 Sl-years per century to determine the unit 
arcsec cy“^. 

The first order post-Schwarzschild correction to 

the perihelion advance can now be calculated from for¬ 
mula ([79]) . For a cosmological constant of A = 10“^^m“^ 
we obtain for the parameters which appear in the expan¬ 
sion ([79]) 


r]i = -.2617993700131308, 
p'‘^{uo) = -1.6972622347915708- 10-^^ 
p"(22o) = 1.3312392184539657- 10“^ 

C{uo) = 1.01210931961465842. 

This leads to a correction of 

^toTT = ^-^2 - 10“^^ arcsec cy“^. (83) 

This result also compares well to [2^ where the perihel 
advance of Merkur does not change within the given ac¬ 
curacy when considered in Schwarzschild-de Sitter space- 
time. The value of the correction is also far beyond the 
measurement accuracy of 0.002 arcsec cy“^ for the peri¬ 
helion advance of Mercury. 

However, for an extreme case the influence of the cos¬ 
mological constant on the periastron advance may be¬ 
come measurable. The orbital data of quasar QJ287 re¬ 
ported in 111,113 indicates that the correction to the pe¬ 
riastron advance will be some orders of magnitude 

larger than the correction in the case of Merkur. Indeed, 
when we calculate from this data the energy parameter 
p and the angular momentum parameter A, 

/i = 0.982166, A = 0.092317, 


(84) 
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we obtain 


~ 10 arcsec cy ^. (85) 

VIII. GEODESICS IN HIGHER DIMENSIONAL 
SCHWARZSCHILD SPACE-TIMES 

We want to show here that our method for solving 
the equation of motion in Schwarzschild-(anti-)de Sit¬ 
ter space-times can also be applied to solve the geodesic 
equation in, e.g., higher dimensional Schwarzschild 
space-times. The metric of a Schwarzschild space-time 
in d dimensions is given by 0 

- 2 , ( 86 ) 

where cM\ = and = dOi + sin^ OidVif for i > 1. 

Because of spherical symmetry, we again restrict the con¬ 
siderations to the equatorial plane by setting Oi = ^ for 
all i. With the conserved energy E and angular momen¬ 
tum L as well as the substitution u = ^ the geodesic 
equation reduces to 

7 \ 2 

^—u‘^— = Pd-i{u), ( 87 ) 

E 9 

where the parameters X = and ji = have the 
same meaning as in the Schwarzschild-(anti-)de Sitter 
case m- For d = 4 this equation reduces of course to 
the Schwarzschild case 0 . For d = 5 a substitution 
u = ^ E I where I is a zero of P4 reduces the differential 
equation (j87|) to 



FIG. 10. Arrangement of zeros of the polynomial Pd-i for 
d = 6. The gray scale code is the same as in Fig. [5] 


where (fe^e = Here cpi is chosen in such a way 

that {2uj)~^(fQ^Q is an element of the theta divisor 0^^ 
and (Pq = ipo dzi depends only on the initial values 

1^0 and cpo. 

The case d = 7 corresponds to a polynomial Pq of 
degree six. If we apply a substitution u = ^ P I where I 
is a zero of Pe, we obtain the differential equation 

+ • • • (91) 

with some constants bi. This can be solved in exactly the 
same way as the differential equation m- The solution 
is 




rs 

u{(p) 


rs—f -r, 

^n^0,7) 


(92) 


ip) ^ + • • • + box° . ( 88 ) 

With an additional substitution x = ^ (4^ ~ ^) 
equation acquires the form (|6TTl which can be solved by 
Weierstrass’ elliptic functions. 

In the case of a d = 6-dimensional Schwarzschild 
space-time, however, the differential equation (j87j) com¬ 
prises a polynomial of degree five on the right-hand side. 
This now can be solved by means of our method. The 
only difference is that the physical angle (p is now given 
by 

r dv' 7^ 

^-^0= / - 7 ^= dz, (89) 

uq y Pbip ) d uq 


what corresponds to dzi rather than to dz 2 as it was in 
the Schwarzschild-(anti-)de Sitter case. This means that 
the solution of the geodesic equation in six-dimensional 
Schwarzschild space-time is given by 


r{ip) = 


rs 

u{ip) 


rs —-r , 


(90) 


where 7 = (^ ^-Po j again, cpi is selected such 

that {2uj)~^ip^ ^ is an element of the theta divisor 0^^ 
and (p'q — ip qP dz 2 . The only difference to the solution 
of the geodesic equation in Schwarzschild-(anti-)de Sitter 
space-time is that the periods uj and uo' and, hence, the 
matrix r will be different due to the different coefficients 
in the polynomials appearing on the right-hand side of 
each differential equation. 

Figure [To] shows the arrangement of zeros for the six¬ 
dimensional Schwarzschild space-time. The gray scale 
code is the same as in the Schwarzschild-(anti-)de Sitter 
case. There are no periodic bound orbits for any values 
of /i and A since the polynomial Pq possesses at most two 
positive zeros. Some resulting orbits for chosen parame¬ 
ters /i and A from different regions in Fig. [10] are shown 
in Fig. [TTl 

With our method it is also possible to analyti¬ 
cally calculate the orbits of particles and light rays in 
Schwarzschild and Schwarzschild-(anti-)de Sitter space- 
times of up to 11 dimensions and in Reissner-Nordstrom- 
(anti-)de Sitter space-times of up to 7 dimensions. Cor¬ 
responding work is in progress |42| . 
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(a) Terminating orbit for = 1.1, A = 0.3 (b) Escape orbit for ij. = 1.1, X = 0.3 (gray 

(gray region in Fig. 1101) . region in Fig. 1101) . 



(c) Terminating orbit for = 0.9, A = 0.3 (d) Infinite terminating orbit for ii = 1.5, 

(light gray region in Fig. HOD . A = 0.6 (white region in Fig. 1101) . 


FIG. 11. Orbits for chosen values of /i and A in six-dimensional Schwarzschild space-time. The black circle indicates the 
Schwarzschild radius. 


IX. SUMMARY AND OUTLOOK 

In this paper the explicit analytical solution for the 
geodesic motion of a point particle in a Schwarzschild- 
(anti-)de Sitter space-time has been presented. We were 
able to explicitly determine all possible solutions and to 
classify them. 

Analytic solutions are the starting point for approxi¬ 
mation methods for the description of real stellar, plan¬ 
etary, comet, asteroid, or satellite trajectories (see e.g. 

Analytic solutions of the geodesic equation can 
also serve as test beds for numerical codes for the dy¬ 
namics of binary systems in the extreme stellar mass ra¬ 
tio case (extreme mass ratio inspirals, EMRIs) and also 
for the calculation of corresponding gravitational wave 
templates. 

The methods presented here are not limited to the 
Schwarzschild-de Sitter case. They also can be applied 
to higher dimensional space-times like Schwarzschild- 
(anti-)de Sitter space-times of up to 11 dimensions and 


Reissner-Nordstom-(anti-)de Sitter space-times of up 
to 7 dimensions. Also for space-times as general as 
Plebahski-Demiahski without acceleration the developed 
method can be applied [i^. In this case polynomials of 
6th order appear which will slightly complicate the struc¬ 
ture of the orbits. It should be noted that this method, 
however, is not capable to solve equations of motions with 
an underlying polynomial of 7th or higher order. In such 
cases we have to enlarge the number of variables to be 
three or more. Then the Abel map is a mapping between 
3 or higher dimensional spaces and it is not clear how to 
constrain this mapping in order to reduce the number of 
variables. 
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Appendix A: Explicit computation of the analytical 
solution 

We describe now in detail the explicit computation of 
the analytical solution derived in Sec. IV. According to 
m the parameters A and p are related the energy 
F, angular momentum L and cosmological constant A 
which gives the polynomial P5 in (p! 3 ]) . The zeros of P5 
have to be determined numerically using, e.g., a Newton 
method. The zeros of P5 already characterize the type of 
orbit as described in Sec. V. 

The Riemannian surface corresponding to is given 
by a two-torus, see Fig. [3l The fundamental paths ai, a2 
and hi , 62 are always the same on this two-torus but 
change in the complex plane depending on the config¬ 
uration of branch cuts, which have to chosen according 
to the zeros of P5. There are eight principally different 
arrangement of zeros of P5 which are shown in Fig. [121 
In each case we choose the branch cuts in such a way 
that we have a maximum number of pure real branch 
cuts. The branch cuts and the resulting canonically fun¬ 
damental paths are also shown in Fig. [ 12 ] together with 
the appropriate sign of the square root on one sheet (with 
reversed sign on the other sheet). 

Next we have to calculate the period matrices (2cc;, 2 uj') 
and (2?7, 2r]') given in (|M|) and, from that, the normalized 
period matrix r defined in (| 27 j) . The periods correspond¬ 
ing to paths Ui, i = 1 , 2 , can be obtained by integration 
from e2i-i to 62^; the integration back from e2i to 62^-1 
yields just the same value due to the different sign of 
the square root. The periods corresponding to paths 
j = 1 , 2 , can be calculated by integrating from 62j to 
62^+1. For paths which encircle one or two complex ze¬ 
ros, the most convenient way to calculate the integral is 
to take a path running straight from the complex zero to 
the real axis and then to proceed along the real axis. Of 
course, fundamental paths can be distorted but only in 
such a way that they do not cross other branch cuts or 
paths as the one showed in Fig. [121 

The branch points are always singularities of the in¬ 
tegrand and, thus, may cause numerical problems. These 
problems can be handled by a partial integration of the 

I 


type 


1: 


dx 


= 2 




- Cj) 


^/Q{x) 


-f 

J a 




VQi^y 


Q'{x)dx , 
(Al) 


where Q{x) = ~ ^3) P some point on the 

integration path. As a consequence, the resulting inte¬ 
grand is no longer singular. In the case that all five roots 
of P5 are real, the computation then is straightforward 
since P5 is always real on the integration paths. If some 
roots are complex, the integration is more involved and 
one has to take into account that the real part of P5 is 
symmetric with respect to the real axis while the imag¬ 
inary part is antisymmetric. One also has to carefully 
choose the appropriate branch of the square root. 

After having calculated all the periods we choose ini¬ 
tial values 1^0 and po and determine = po ^ dz2- 
Note that the integral is just a sum of the half-periods 
already calculated above provided uq, as usual, 
is taken to be one of the zeros of P5. For every (p we 
have to find now the dummy parameter pi such that 
{2uj)~^(fQ is an element of the theta divisor 0^ = 

{z|'d (^1/2)5(172) This dummy param¬ 

eter has no physical meaning and depends, beside (7?, 
on the initial values uq, (po and the normalized pe¬ 
riod matrix r. The value of (pi can be computed with 
a Newton method for the function g o h : (D ^ (D, 

if A A 'd At). For 

p = (po and 1^0 being a branch point, the dummy param¬ 
eter (pi is explicitly known. From (p — (p^ = dz2 and 

P 


ro/2\ 

(0 

[A/y 

. U/2 j 


obtain in this case pi = dzi. 

After these calculations we are finally able to compute 


/ N rs 0-2 Ae) 

= — 7 -^ = -rs- 


uA) 


(JiAe) 


(A2) 


The solution r{p) can be computed pointwise for any 
value of p with, in principle, arbitrary accuracy. 


Appendix B: Calculation of post—Schwarzschild 
period 

Toward an integration of the second term of (HH) we 
first obtain 


/ 

J ai 


-h A 


-dx = 


(4z+i)2 + A 


a, x^P^{x)PpAA) Ja (4z + i)2. 4^Pw{z)A4^Pw{z) 


Mz 


1 

42 


dz 


^ ^ _ 

Ja {4z + 1)^Pw{z)PpA(z) 


A Pw{z)\/Pw{z) 


(Bl) 
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zeros (ei, 02) 




04) 


FIG. 12. Branch cuts (fat lines) and fundamental paths in the complex plane for all arrangements of zeros ei,...,e5 of a 
polynomial of degree 5 and ee = 00. The branch cuts are drawn from e2z-i to e2z. The completion of the 5 -paths on the other 
sheet is indicated by dashed lines. 


In the following we will represent the functions 

1 _ . . 1 


Fi{z) = 


Pw{z) ’ 


F2{z) = 


{Az + \YPw{z) 


(B2) 


as linear combinations of the Weierstrass elliptic function 
p as well as the Weierstrass C function. The reason is that 
these functions can be integrated easily since C' — ~P 
and (logcr)' = C, where a is the Weierstrass’ a-function: 


f p{u - uo)du = Ci'fiO) - uo) - Cilii) - uo) (B3) 
/ C{u - uo)du = log (7(7(1) - Uo) - logcr(7(0) - Mo) , 

J7 


(B4) 


where 7(0) = ui and 7(1) = ui uji as above (the 
branches of log will be discussed later). For the and 


cr-functions we have a quasiperiodicity 


C{u + U}j) = C{u)+Vj 

a(u + ujj)= e’'i(“+“^/2)+^V(M), 


where r]j are periods of second kind given by 


Vi = - 

m = - 


i 

i 


zdz 

\J Pw (z) 

zdz 

\/Pw{z) 



p{u)du 


p{u)du. 


(B5) 

(B6) 


(B7) 
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Thus, JEHl) and (lB4l) can be rewritten as 


In summary, we obtain 


/ p{u-uo)d 
/ C(m - Mo)<^ 


u = -r]i 

u = 771 ( 1^1 — uq \ooi) + Tri + 2 iTik 


(B 8 ) 

(B9) 


with k ^ Wj. 


1. Integration of Fi/\/Pw• 


First we substitute 2: = p{u) 

1 1 


Fi{z) 


Pw{z) Pw{p{u)) p'{uy 


=: /iH- 


(BIO) 

The function fi only possesses poles of second order in 
pi, p 2 , and ps. In a neighborhood of pj, the function fi 
can be expanded as 


/l(w) = 


aj2 


aji 


{u - Pjf 


pj 


■ holomorphic part. (Bll) 


Since p'{pj -\- z)‘^ = p'{pj — z)^ for all j and /i is sym¬ 
metric with respect to all pj and, therefore, depends only 
on even powers of {u — pj) so that aji = 0. The constant 
aj 2 can be evaluated with a comparison of coefficients. 
For this, we note that p'{pj) = 0 = p"'{pj) and, thus. 


00 

pju) = p”{pj){u - Pj) + - Pjf (B12) 

2 = 3 


in a neighborhood of pj and for some constants c^. If 
we square both sides of the equation, we see that 
contains only even powers of {u — pj) larger than 1 . It 
follows 


. ,_ — / fi{u)du 

A Pw{z)^JPw{z) 

= / 'paj2{p{u -pj)- Zj)du 

j=i 


J = -L 
3 




(B17) 


2. Integration of F^j^jPw'- 


Again, we first substitute z = p{u) and obtain 


F2{z) = 


1 


(4z + ^)^Pw(z) 

1 


((4p(u) + Pp'(u)) 


(4p(u) + ^)^Pw(p(u)) 

- =; f2(u). (B18) 


The function /2 possesses poles of second order in pi, p 2 , 
P 3 and in all uq e R such that p{uo) = Since we 

assumed that the considered orbit is bound, all zeros of 
P 3 have to be positive and, thus, zi > Z 2 > z^ > 

This means that 0 < Im(7Xo) < Im(p 3 ). The function p 
is even and, hence, also uq \= uj 2 — uq G P is a pole of 
second order (see Fig. [9j). 

Since p is symmetric with respect to pj, the function 
/2 can be expanded in the same way as above as 


1 = fi{u)p'{u)‘^ = aj 2 p"{pj)‘^FhigheT powers of (u-pj ). 

(B13) 

From that it follows aj 2 = T The function 

p{u — pj) has only one pole of second order in pj with 
zero residue. Therefore, the difference 

3 

h{u) - aj 2 p{u - Pj) (B14) 

i=i 

is a holomorphic elliptic function and, thus, is constant. 
This yields 

3 

hiu) = Y “i2p(w - pj) +Cl. (B15) 

7 = 1 

The constant ci can be determined by /i(0) = 0 using 
the relation p{—pj) = p{pj) = Zj: 

3 

ci = -^a^- 2 ^i. (B16) 

7=1 


f 2 {u) = 7 -T holomorphic part (B19) 

[u — Pj ) 

in a neighborhood of pj. An expansion of {4:p{u) + 
^)p\u) near pj yields 

{ 4 p{u) + i) p'{u) = aji(u - pj) + aj2{u - pj)"^ 

+ higher order terms (B20) 

because of p'{pj) = 0. The coefficients are given by 

= {{Mu) + l) p'iu))Y^ 

= {4zj + i) p'jpj) (B21) 

a ,-2 = ((4p(u) + i) p'{u))Y^ = 0 . (B22) 

A comparison of coefficients 

l = /2(u)((4p(u) + i)p'(u))" 

= aj 2 a‘ji + higher powers of {u — pj) (B23) 
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yields 


In a neighborhood of fto, the function /2 is given by 


+ holomorphic part. (B24) 

in a neighborhood of pj. 

The same procedure will be carried through for uq and 
uq. We have 

f 2 {u) = 7—^ H- — -h holomorphic part (B25) 

[u — Uq)^ u — Uq 


f 2 (u) = 7—H-+ holomorphic part. (B32) 

[U — UqY U — Uq 

As p' is an odd and p" an even function we get for the 
coefficients of the expansion of {^p{u) + ^)p'{u) near uq 
with Uq = UJ 2 — Uq the relations 

^i=/3i, h = -P2 (B33) 

and, therefore, 

6i = -6i, b2 = b2. (B34) 


and 


(4p(it) + i) p'{u) = Pi{u - Mo) + hiu - Uof 

+ higher order terms (B26) 

near uq. The coefficients of (jB26p read 

/3i = 4 p'(mo)2 (B27) 

P2 = 6p'(mo)p"(mo) . (B28) 

Again, a comparison of coefficients 

1 = / 2 (w) ((4 p(m) + I) p'{u)f 

= hPl + (262/3 i /32 + hiPl){u - Uq) 

+ higher order terms (B29) 

yields 

16p'(mo)^ 

61 = -2f3,f32b2f3p = 

16 p'(mo)5 


(B30) 

(B31) 


Summarized, the function 
3 

92{u) :='^aj2p{u - pj) + b2{p{u - uq) + p(m - mq )) 
i=i 

+ bi{C{u - Uq) - C{U - Uq)) (B35) 

has the same poles with the same coefficients as /2. 
Therefore, /2 — ^2 is a holomorphic elliptic function and, 
thus, is equal to a constant C2. This constant can be de¬ 
termined by the condition 0 = /2(0) = ^2(0) + C2 which 
yields 


3 

C2 = - ^aj2Zj + \b2 - bi{C{uQ) - C{uo)) ■ (B36) 

j=i 

As a consequence, 

3 

f2{u) = y] aj2ip{u - Pj) - Zj) 

i=i 

+ b 2 {p{u - Uq) + p(m - Mo) + g) 

+ bl{C{u - Uq) - C(m - Mo) + C{uo) - ({uq)) . 

(B37) 

Now we can carry through the integration of the second 
term on the right-hand side of m- 
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/ ^2(2:)^=== = / h{u)du 

■I A ^JPw{z) Jj 

y~^ aj2 p(m - pj)du - ZjUi^ y ~ '*’ ~ wo)'^w + 5^1 

bi / C(m - Mo) - C(^i - wo)'^w + &iWi(C(wo) - C(wo)) 

J-i 


J = ^ 
+ 1 


= - + ^j^i) + h (1^1 - 27/i) + 6i(r/i('Uo - Uo) + 27ri(/ci - /C 2 )) 

i=i 

+ 6iu;i(C(^o) -C(^o)) 

3 

= — CLj2{r]i + + 62 (|<^i ~ 2771) 

i=i 

+ 6i(7/i(u;2 - 27Xo) + 27ri{ki - fe)) + biuJi{2C{uo) - 772 ) • 


(B38) 


The difference {ki — k 2 ) can be calculated as follows. 
First note that via x = ^p{u) ^ uq corresponds to 
0 and 7Xi to X3. Since 0 < X3 for bound orbits under 
consideration we have Ini(77o) < Ini(7Xi) < Ini('Lto). Let 
now I be determined by 

/ C{u - uo)du = r]i{u2 - fto + —) + 7ri + 27ril , 

Ju2 ^ 

(B39) 

where U 2 G TR is such that Im(7X2) > Im('Uo) > Im(7Xi) > 
Im(7Xo). From 


(i) I does not depend on uo as long as Im('Lto) < Im(7X2) 


holds and, thus, 

CJi 

/ C,{u—uo)du = irii{u2—U{)^-—)^T:i^2T:il (B40) 

J U2 ^ 

and 

(ii) (jBdOp holds also for U 2 replaced by ui by Cauchy’s 
integral formula for the rectangle with corners 7x1, 
7Xi + cji, 7x2 + uJi and 7x2. 

it follows that I = ki. 

We show now that k 2 = / + 1 and, thus, k 2 = ki -\- 
1. Consider the counterclockwise oriented rectangle with 
corners 7x1 , txi+cji, 7X2+Cc;i and 7x2. Let c be the boundary 
of this rectangle but with a two symmetric small bumps 
such that c encircles the pole txq of C('^ “ ^ 0 ) with residue 
1, but not TXo + cji. Then the residue theorem gives 


27rx = 


^ ({u — uo)du 

nUi-\-UJi pU2-\-UJi pU2 nU\ 

= / ("(tX — 7Xo)d7X + / ("(tX — 7Xo)d7X + / C(7X — 7Xo)<^'^ + / (^{u — UQ)dVi 

J Ui J U2-\-UJi J U2 

nU\+UJ\ PU 2 +UJ 1 

= / ({u-uo)du- / ({u - uo)dur]i{u 2 - Ui) 

J J Uo 


= - fio + \^i) + ttx + 27 Tik 2 - (771(7x2 - TXo + \^i) + ttx + 2 ml) + 771(7x2 - txi) 

= 2 m{k2 — 1) . 

With the Legendre relation 77iCo’2 — ^2<^i = 27ri we finally obtain 


£f2 (^) 


dz 


771 + ZjLUi 


\/Pw{z) ^ \ (Azj + p"{pjY, 


1^1 - 2t?i _ 6 p"{uo) 
16 p'(mo)^ 16 p'(uo)® 


{-piuo + a;iC(Mo)) • 


(B41) 


(B42) 


Note that though the values p'(7Xo)^, txq and C('^o) appearing in the last part of the right-hand side are all purely 
imaginary, the hole term is real. 
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